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Abstract. The spectral properties of Dirac operators on (0, 1) with potentials that 
belong entrywise to L p (0,l), for some p € [1, oo), are studied. The algorithm of 
reconstruction of the potential from two spectra or from one spectrum and the cor- 
responding norming constants is established, and a complete solution of the inverse 
spectral problem is provided. 



1. Introduction 



The main aim of the present article is to solve the direct and inverse spectral problems 
for one-dimensional Dirac operators on a finite interval under possibly least restrictive 
assumptions on their potentials. Namely, the Dirac operators under consideration are 
generated by the differential expressions 



dx 

and some boundary conditions, where 

an b- ( l \ o(x)- ( qi ^ x) q2 ^ x) 



and qi and q 2 are real-valued functions from L p (0, 1), p G [1, oo). To simplify unessen- 
tial technicalities, we shall only consider the boundary conditions that correspond to 
the Neumann-Dirichlet and Neumann ones in the case of Sturm-Liouville equations, 
although other boundary conditions can be treated in a similar manner (cf. the study 
of Sturm-Liouville operators with nonsmooth potentials and various boundary con- 



ditions in [261 152]). The corresponding Dirac operators srf\ and stf^ in the Hilbert 
space H := L 2 (0, 1) x L 2 (0, 1) act according to the formula = £qU on the domains 

dorn^- := {u = (u 1 ,u 2 f | u l7 u 2 G AC(0, l),£ Q u G M,u 2 (0) = %(1) = 0}. 

It is well known that the operators sf\ and srf 2 are selfadjoint in H and have 
simple discrete spectra accumulating at — oo and +oo. Our primary goal is two-fold: 
firstly, to give a complete description of the spectra of srf\ and s$ 2 for potentials Q of the 
form (11. ip with qi,q 2 G L p (0, 1) for some p G [1, oo) — i.e., to solve the direct spectral 
problem, — and, secondly, to give an algorithm of reconstruction of these operators from 
their spectra or from one spectrum and the corresponding norming constants — i.e., to 
solve the inverse spectral problem. 
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Ever since P. Dirac suggested in 1929 the equation (later named after him) modelling 
the evolution of spin-| particles in the relativistic quantum mechanics [51], its range 
of applicability in various areas of physics and mathematics has been continuously 
expanding. In particular, in 1973 Ablowitz, Kaup, Newell, and Segur [1] discovered that 
the Dirac equation is related to a nonlinear wave equation (the "modified Korteweg-de 
Vries equation", a member of the AKNS-ZS hierarchy, see [2j [56]) in the same manner 
as the Schrodinger equation is related to the KdV equations, and this stimulated the 
increasing interest in direct and inverse problems for Dirac operators in both physical 
and mathematical literature. Earlier in 1966, Gasymov and Levitan solved the inverse 
problems for Dirac operators on M + by using the spectral function [16] and by the 
scattering phase [15]. Their investigations were continued and further developed in 
many directions. The reference list is so vast, that we can only mention those papers, 
which, in our opinion, are most pertinent to our topic and refer the reader to the 
bibliography cited therein for further material. The books by Levitan and Sargsjan [36] 
and by Thaller [31] may serve as a good introduction to the (respectively mathematical 
and physical part of the) theory of Dirac operators. 

The inverse scattering theory was developed for Dirac operators on the axis in [TT) 
[T2"| |2~T| |2"3"1 EES] , for Dirac systems of order 2n on semiaxis in [T3] , and for more gen- 
eral canonical systems on IR in [50J. The nonself adjoint case was treated in [37] and 
nonstationary scattering, including point interactions, in [17] and [1] , respectively. Re- 
construction from the spectral function on semiaxis was done in [51] for a general 
boundary condition at x = and in [53] in the case of an interface condition in an 
interior point; the general first order systems in L2(R+,C 2n ) were recently treated 
in [33]. Inverse problems in the periodic case were studied in [291 EI], an d the Weyl- 
Titschmarsh m-function was used to recover the potential of the Dirac operator in [5] 
and of the Dirac systems of order 2n in [7J [T7J 09] (see also the detailed reference lists 
therein) . 

The inverse problems for Dirac operators on a finite interval have also been studied 
in detail. Reconstruction of a continuous potential from two spectra was carried out 
in [H], from one spectrum and the norming constants (in the presence of a Coulomb- 
type singularity) in [9], and from the spectral function in [ID]. Explicit formulae for 
solutions (based on the degenerate Gelfand-Levitan-Marchenko equation) in the case 
where finitely many spectral data are perturbed were given in [8]. Uniqueness results 
for other types of inverse problems were established — e.g., for mixed spectral [25] or 
interior [12] data, nonseparated boundary conditions [13], or for the weighted Dirac 
equations [55]. Ambartsumyan-type theorems were proved in [21] and for the matrix 
case in J3UJ- Finally, uniqueness of the inverse problem for general Dirac-type systems 
of order 2n was recently established in [3H1 13"9"] . 

We observe that in the above-cited papers the inverse spectral problems for Dirac 
operators on a finite interval were considered for continuous potentials only, which 
excludes, e.g., the important case of piecewise constant potentials. We remove this 
restriction by allowing potentials belonging entrywise to L p (0, 1), p G [l,oo), and 
completely solve the inverse spectral problem for Dirac operators in this class (see 
Theorems 13.11 13.31 15.11 and 15.71) . The main idea of the proof rests on the fact that 
the transformation operators for Dirac operators under consideration satisfy not only 
the classical Gelfand-Levitan-Marchenko equation (14.41) . but also its counterpart (14.91) . 
which was used by Krein [32] in the study of the inverse problem for impedance Sturm- 
Liouville equations (see also [3], lOj). This "Krein equation" survives the passage to 
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the limit in the L p -topology and thus allows us to treat potentials belonging to L p (0, 1) 
entrywise. 

The paper is organized as follows. In Section [2] transformation operators are con- 
structed and some of their properties are established. Based on this, in Section [3] we 
find the asymptotics of eigenvalues and norming constants for the operators and 
&/ 2 - The Gelfand-Levitan-Marchenko and Krein equations, which relate the spectral 
data and the transformation operators, are derived in Section HJ and the solution of the 
inverse spectral problem is given in Section [5j Finally, two appendices contain some 
facts related to harmonic analysis and the factorisation theory in operator algebras. 

Throughout the paper, we shall denote by (•, •) the scalar product in HI and by <M 2 
the algebra of 2 x 2 matrices with complex entries endowed with the operator norm 
| • | of the Euclidean space C 2 . Where no confusion arises, we abbreviate L p (0, 1) to L p 
and write L p (^f 2 ) for the space L p ((0, l),^#2j of ^2-valued functions on (0, 1) with 
complex-valued entries and the norm 

\\V\\ Lp := (j\v{t)\vdt) llP . 

Also, (x, yY shall stand for the column-vector in C 2 with components x and y. 



2. Transformation operators 

Assume that Q G L p (^ 2 ), P G [l,oo), is of the form (11 .11) and denote by U(-) = 
U(-,X) the Cauchy matrix corresponding to the equation £qu = Au. In other words, 
U is a 2 x 2 matrix-valued function satisfying the equation 

(2.1) B < ^ + QU = XU 

dx 

and the initial condition U(0) — I := diag(l, 1). Denoting by c(-, A) := (ci(-, A), c 2 (-, A)) 
and s(-, A) := (si(-, A), s 2 (-, A))* the solutions of the equation £qu = Au satisfying the 
initial conditions ci(0, A) = ^2(0, A) = 1 and 02(0, A) = si(0, A) = 0, we find that 

U(x,\) = N X ^\ S f^\ 

\C2{x, A) s 2 [x, A) 

Our next aim is to derive an integral representation for U of a special form. 
Theorem 2.1. Assume that Q 6 L p (^2), P £ [1, 00). Then 

(2.2) U{x, A) = e~ XxB + [ e- x{x - 2s)B P(x, s) ds, 

Jo 

where the matrix-valued function P = Pq has the following properties: 

(a) for every x G [0, 1] the function P(x, ■ ) belongs to L p (^ 2 ); 

(b) the mapping Vq : x 1— > Pq(x, ■ ) G L p (^ 2 ) is continuous on [0, 1]; 

(c) the function Vq depends continuously in C([0, 1], L p (^ 2 )) on Q G L p {^ 2 ). 

Proof. The standard variation of constant arguments show that U satisfies the equiv- 
alent integral equation (recall that B 2 = —I) 

(2.3) U{x, A) = e~ XxB + / e- x{x - t)B BQ(t)U(t) dt, 

Jo 
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which can be solved by the method of successive approximations. Namely, with 

(2.4) U Q (x) := e~ XxB and U n+1 (x) := f e~ x{x - t)B BQ(t)U n (t) dt for n > 0, 

Jo 

the solution of (12.31) formally equals J2^=o U n - Assume that we have proved that 

oo 

(2-5) ^2\\U n \ 



loo < OO, 



n=0 



where ||£/ra||oo := sup^^y |f/ n (a;)|. Differentiating then the recurrence relations (12.41) . 
we find that 

U' n {x) = -XBU n (x)+BQ(x)U n ^(x), 

which in view of (12.51) shows that the series ^^1 U n converges in the topology of 
the space ((0, 1), ^ 2 ) to some ^ 2 -valued function V. This function V solves (12. ip 
and satisfies the initial condition V(0) = /, and hence it coincides with the Cauchy 
matrix U . 

To justify (12.51) . we use the identity 

(2.6) e~ XxB Q(t) = Q(t)e XxB , x, t e [0, 1], 
in the recurrence relations (12 .4p and derive the formula 

(2.7) U n (x) = [ e -^-2Mt))B BQ{fi) . . . BQ{fn) ^ dtni 

Jn n (x) 

in which we have set 

II n (x) = {t := (t u ...,t n )eW n \0<t n <---<ti<x}, 

n 

en(t) = D-iy+v 
i=i 

Upon the change of variables s = £ n (t), yi = tz+i, I = 1,2, ... ,n — 1, we recast the 
integral in (12. 7p as 

U n {x) = [ X e~ x ^ B P n (x,s)ds, 
Jo 

where P±(x, s) = BQ(s) and, for all n G N and < s < x < 1, 

(2.8) P n+1 (x, s)= [ BQ(s + Uy)) BQ(yi) ■ ■ ■ BQ(y n ) d Vl ... dy n , 
with 

U* n (x, s) = {y = (yi, ...,y n )eR n \0<y n < y n ^ <---< yi <s + £ n (y) < 4- 

For convenience, we extend the functions P n , n > 2, to the whole square [0, 1] x [0, 1] 
by setting P n (x, s) — for < x < s < 1. 
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Using the Holder inequality and Fubini's theorem, we find that, for every n G N, the 
function P n+1 (x, ■) belongs to L p (^ 2 ) and that 

\\Pn + l(x, -)\\ P Lp = [ \P n+1 (x, S )\^ds 

Jo 

<{n\) 1 ^ [ [ \Q(s + ^- 1 (y))\ p \Q(yi)\ P ---\Q(yn)\ p dy 1 ...dy n ds 
= {niy- p I \Q(t 1 )\ p ---\Q(t n+1 )\Pdt 1 ...dt n+1 

/ f x \ n+1 

(/ior) < 



i f r p y +1 \\Q\\ ( C 1)P 

(n\)p{n + 1) \J |g V - (n!)f ' 



Henceforth with C := max^r-i^i |e XxB \ we have 

WQWl 



U n (x)<C P n (x,s)ds<CP n (x,-)\ Lp <C, 

Jo (n-1)! 

and (12.51) follows. 

Moreover, the above inequality implies that the series Y^Li P n ( x , •) converges in 
L p {^2) to some function P(x, •) and yields the estimate 

(2-9) ||P(x, -)lk <E^W = HOlk«P{llOHU 

n=l ' 

for all x G [0, 1]. This establishes (a). 

Assume that Q is another potential in L p (^ 2 ) and denote by P n the corresponding 
functions constructed as above but for Q instead of Q; then similar calculations on 
account of the inequality 

n n n ^ 

| JJctfc - Y[b k \ P < K - h\ JJ(|flil + \bj\)j 

k=l k=l k=l j^k 

n 

k=l j^k 

lead to the estimate 

/ in _|_ 1 \ p / _\np 

(2.10) ||p t+1 (x,.)-P« + i(x,-)lk< (-±-) + • 

It follows that 

(2.11) ||P Q (x, • ) - Pq(x, ■ )\\ Lp < (1 + 2r)e 2 mg - Qlk 

as soon as r is such that ||<3||z, , ||Q||l p < r - 

Observe that if Q G C([0, 1],^ 2 ), then the functions P n , n > 2, are continuous in 
the square [0, 1] x [0, 1], and, moreover, 

max \PJx,s)\ < 



nV""j ° / — i ) 



6 S. ALBE VERIO, R. HRYNIV, AND YA. MYKYTYUK 

so that the function [0, 1]9ih Pq( x , ■ ) £ L p (^ 2 ) is continuous. Since the potential 
Q G L p (^ 2 ) is the limit in L p (^# 2 ) of potentials Q n G C([0, 1], -#2), estimate (12. lip 
yields both assertions (b) and (c). The proof is complete. □ 

Corollary 2.2. Assume that Q e L p (^# 2 ) one? se£ 

00 00 



2n-lj 

n=l n=l 



where the functions P n are given by formula (12.81) . Set c (x, A) := (cos Ax, sin Ax)* and 

P(x,t) = Pg(x,t) := P + (x,t) + p-(x,t)J, 
{2AJ) K{x, t) = K Q (x, t) := \ [R(x, 2^) + P(x, 2±*) J] , 

where J = diag{l, —1}. Then the vector-function c( • , A) is given by 

(2.13) c(x, A) = cq(x, A) + / K(x,t)c (t,X)dt. 



Proof. Using (12.61) and (12.81) . we conclude that 

e~ x ^ B P 2n (x,s) = P 2n (x,s)e~^^ B , 

Therefore equality (12.21) can be written as 

U(x, A) = e- XxB + /" P + (x, s)e- A(x - 2s)B ds + [* P~(x, s)e x{x ~ 2s)B ds. 
Jo Jo 

Observing that 



-XxB _ I cos Ax — sin Xx 
sin Ax cos Ax 



and taking the first column of the above equality, we get 

PX 

c(x, A) = c (x, A) + / R(x, s)c (x — 2s, A) ds. 
Jo 



Since 



x/2 ! 

2 J 



R(x, s)c (x — 2s, A) ds = - / R(x, ^)c (t, A) dt 



/•X fX 

\ R[x, s)c (x — 2s, A) ds = / R(x, s) Jc (2s — x, A) ds 

Jx/2 Jx/2 



c/2 Jx/2 

-x 

R(x,^)Jc (t,\)dt, 



1 

~ 2 

the required relation follows. □ 
Equality (12.131) shows that the operator J? + defined by 



fX 

{J + JT)u(x) = u(x) + / K(x, t)u(t) d£ 



transforms the solution of the equation l^u = Au (i.e., with a potential Q equal to 
zero identically) subject to the initial conditions «i(0) = 1, u 2 (0) = into the solution 
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of the equation £qu = Au satisfying the same initial conditions. Denote by s^q the 
operator in HI acting as Aqw = £qu on the domain 

domK?) := {u = (ui, u 2 y G W^(0, 1) x W 2 \0, 1) | u 2 (0) = 0}; 

then + is in fact the transformation operator for £$q and j&q, i.e., £^q{^ + Jff) = 
(J? + Jf)gf , see Theorem E3 

The operator JgT possesses some important properties, which we now establish. De- 
note by G p {^2) the set of measurable 2x2 matrix-valued functions K on [0, 1] x [0, 1] 
having the property that, for each x and t in [0, 1], the matrix- valued functions K(x, ■) 
and K(-,t) belong to L p (^ 2 ) and, moreover, the mappings 



[0, 1]3ih K(x, •) G L p (^ 2 ), [0, 1] 3 t i-> K(-, t) G L p 



2) 



are continuous (i.e., they coincide a.e. with some continuous mappings from [0, 1] into 
L p {^ 2 )). The set G p (^#2) becomes a Banach space under the norm 

(2.14) \\K\\ Gp :=max{ max \\K(x, -)\\ Lp , max :\\K(;t)\\L p }. 

a;G[0,l] te[0,l] 

We also denote by W v {^ 2 ) the set of the integral operators in EI with kernels K 
from G p (^ 2 ). Under the induced norm \\J^\\& P '■— \\K\\g p , the set C S P {^ 2 ) becomes 
an algebra. The algebra ^ p {^ 2 ) is continuously embedded into the algebra 3§{M) of 
all bounded operators in HI since the functions K belonging to G p {^ 2 ) have finite 
Holmgren norm [22J; moreover, for J?T G C S P [^.M 2 ) the inequality \\J(f\\gg(wi) < \\^\\% 
holds true. 

Theorem 2.3. Assume that Q G L p (^ 2 ); then the integral operator J(f = Jtq with 
kernel K of f)2.12p belongs to @ p (^ 2 ) and, moreover, the mapping L p (^ 2 ) 3 Q l— > 
G ^,(^#2) continuous. 

Proof. In view of relations f)2.12p . we have 

(2.15) K(x, t) = ^ [P+(x, 2^) + P-(x, ^) J + P+(x, ^) J + 2±*)] , 

and hence it suffices to prove the assertions of the theorem for the operators with 
kernels P + (x, ^^) and P~(x, ^). Since the proof is analogous for all four functions, 

we shall give it for the function P + (x, ^) =: P(x,t) only. 

It follows from the proof of Theorem 12.11 that the function P + enjoys the properties 
(a)-(c) of that theorem. Simple arguments based on the change of variables justify the 
validity of the properties (a)-(c) for the kernel P. It thus remains to establish similar 
properties of P with respect to the variable t. 

Assume first that Q is continuous. Changing the variables r\ = s + £(y), y\ = y 2 , 
. . . , y n -i = y n in integral (I2.8p . we arrive at the relation 

BQ(v)P n (v,V - s ) d V, 

which yields 

P+(x,s)= [ BQ(r))P-(r],r}-s)dr). 
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It follows that 



p 

dx 



\\H; t) \\l p = J | J x ^ BQ( V )P-( Vi V -^)dv 

< J d7i\Q(7i)\*J \P-(r,,r]-^)\ p dx 

< \\Q\\l p max \\P-(v,W Lp < \\Q\\ 2p ^p{p\\Q\\l p }, 

cf. (12 .9p . If Q is another continuous potential, then we find analogously that 
\\P Q (;t) - pQ(;t)\\ P Lp < ^WQ ~ QWl max \\P Q ( V ,-)\\l p 

+ ^\Q\\l pS -jP Q (Vr)-P Q (v,-)\\ P Lp . 

Recalling inequality (I2.10p . we conclude that the function 1 1— > Pq(-, t), which belongs to 
C([0, 1], L p (^ 2 )) if Q is continuous, depends therein continuously on Q G C([0, 1], ^#2) 
with respect to the topology of L p (~-# 2 ). Since the space C([0, 1], ^#2) is dense in 
L p ( t y# 2 ), we show by continuity that, for every Q G L p (^ 2 ), the function Pq(-,t) 
belongs to L p (^f 2 ) f° r every fixed t G [0,1], that the mapping t 1— > Pq(-,t) is con- 
tinuous, and that the continuous L p (.^2) -valued function of t, t 1 — > Pg(-,t), depends 
continuously in C([0, l],L p (^# 2 )) on Q G L p {^ 2 ). 

Summing up, we have shown that the function P = Pq belongs to G p (^ 2 ) and 
depends in G p {^ 2 ) continuously on Q G L p {^ 2 ). This establishes the theorem. □ 

Theorem 2.4. Assume that Q G L p {^ 2 ) and let J{f be an integral operator with ker- 
nel K of (I2.12p . Then J? + X is the transformation operator for the pair stfq and s^q, 
i.e., J Q {J? + Jff) = (j? + Jtf)J . 



Proof. Since X belongs to ^ p {^ 2 ) by Theorem 12.31 and its kernel K is lower-diagonal, 
it follows that J(f is a Volterra operator in H and hence J" + X is a homeomorphism 
ofM. 

Write J Q := {J + X)' 1 ^^ + X). In view of fl2TT3|) . c ( ■ , A) is an eigenvector 
of the operator s^q corresponding to the eigenvalue A for every A G C. Denote by Jzf 
the linear hull of the system {co( ■ , A) | A G C}; then the restrictions of the operators 
&/q and £/q onto Jzf coincide. Since Jzf is a core of (see a similar result in [3l 
Theorem 3.3]) and &/q is closed, it follows that s^q C s$q. It remains to observe that 
s^q cannot be a proper extension of s^q since otherwise s^q — and s^q by similarity — 
would have a two-dimensional nullspace, which would contradict the uniqueness of 
solutions to the equation £qu = Au. Thus M) = £^q, and the proof is complete. □ 

3. Direct spectral problem 

The aim of this section is to perform the direct spectral analysis for the Dirac op- 
erators and stf 2 . The main tool of our investigations will be the transformation 
operators constructed in the previous section. 

Theorem 3.1. Assume that Q G L p (^ 2 ); then the eigenvalues (X n ) n & an d (/-tn)nez of 
s^i and stf 2 respectively can be enumerated so that they satisfy the interlacing condition 

(3.1) A n _i < fi n < A n , n G Z, 
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and the asymptotics 
(3.2) 



A n = 7i(n+ §) +e n (gt), 
fi n = im + e n (g 2 ), 



where gi, g 2 G L p and e n (g) := e 27rmx g(x) dx are the Fourier coefficients of a func- 
tion g. 

Proof. The fact that the spectra of and srf 2 interlace (i.e., that between two consec- 
utive eigenvalues of one operator, there is exactly one eigenvalue of the other operator) 
is well known (see, e.g., p3|). It is easily seen that if we have an enumeration of A n 
and fi n obeying (13. 2p for some g±, g 2 G L p , then (13.11) holds for all n with sufficiently 
large \n\. Since the two spectra interlace, we can permute a finite number of indices 
if necessary in such a way that (13.11) becomes valid for all n G Z. This reordering 
amounts to adding trigonometric polynomials of finite degree to the functions g\ and 
g 2 , and the modified functions g\ and g 2 will remain in L p . Therefore it suffices to 
establish Q . 

The equation £qu = Au subject to the initial conditions /i(0) = 1, f 2 (0) = has the 
solution 

fX 

c(x, A) = c (x, A) + / K{x, s)c (s, A) ds, 



where K =: (kji)? l=1 is the kernel of the transformation operator constructed in Corol- 
lary [272J The numbers A n are zeros of the function Ci(l, A), which, after simple trans- 
formations, takes the form 



(3.3) 



where 



ci(l, A) = cos A + / [ftii(s) cos(As) + ki 2 (s) sin (Ax)] ds 
Jo 



cosA+ / /i(s)e tAa ds 



i 




l,s) -ih 2 (l, s)\, s > 0, 

1, -s) +ik 12 (l, -s)], s<0 



m := 

is a function in L p (—1, 1). 

Analogously the eigenvalues of s$ 2 are zeros of the entire function c 2 {\, A), which has 
the form 

(3.4) c 2 (l,A) =sinA + y f 2 {s)e iXs ds 

for some f 2 G L p (—1, 1). 

The required asymptotics for zeros of ci(l, A) and c 2 (l, A) follows now from [27]. □ 

Definition 3.2. We denote by SD p the set of all pairs {(A n ) ng ^, (/v)nez}) i n which 
(A n ) and (/x n ) are sequences of real numbers that obey the interlacing condition (13. ip 
and the asymptotics (13.21) . 

Theorem 13.11 shows that, for any real- valued Q G L p (^ 2 ), the spectra of the opera- 
tors srf\ and srf 2 form an element of SD P . In the reverse direction, Theorem 15.11 claims 
that any element of SD p is composed of the spectra of Dirac operators s$\ and sd 2 for 
some Q G L p (^f 2 ). The reconstruction algorithm uses in fact the spectrum (A n ) of 
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stf\ and the sequence of corresponding norming constants a n := ||c(-,A n )|| 2 , whose 
properties we are going to study next. 



1-2 



Theorem 3.3. Assume that Q G L p (^# 2 )/ then the norming constants a n = ||c( • , X Tl/ m 
have the asymptotics 

(3.5) a n = l + e n (g), 

where g G L p . 

Set 0(A) := Ci(l,A) and ^(A) := c 2 (l,A). We observe that the functions 0(A) and 
■0(A) can be reconstructed from their zeros as follows 



Proposition 3.4. The following equalities hold: 



(3.6) 0(A) = V.p. J] 



A ra — A 



7r(n+|)' 

n=— oo v ^ ' 



^(A) = (A-/i )V.p. ] [ 



/in - A 



7m 



(where the prime means that the factor corresponding to the index n = should be 
omitted); moreover, the products converge uniformly on compact sets. 

It turns out that and -0 (i.e., that the two spectra (A n ) n ez and (/i n )nez) determine 
the norming constants a n as follows. 

Lemma 3.5. The norming constants a n satisfy the following relation: 
(3.7) a n -■ 



0(A n )^(A n ) 



Proof. Since the system of eigenvectors (c n ) neZ , c n := c(-,A n ), of the operator is 
an orthogonal basis of H, we have 



A)-'u= £ 



A„ — A 



In particular, the residue of this expression at A = A n equals — a n (u, c n )c n . 
On the other hand, (s^i — A) -1 can be calculated as 

r r 1 



(u(t),8+(t,X)) C 2dt 



c(x, A) 



+ 



(u(t),c(t,X)) C 2dt 



J W(X) 
s+(x,A) 



W(X) 



where s + (-, A) is a solution of the equation £q(u) = A 2 u subject to the terminal condi- 
tions Mi(l) = and «2(1) = 1, and W(X) is the Wronskian of the solutions c(-, A) and 
s+(-, A), i.e., 

W(X) := ci(x, X)$2 (x, A) - c 2 (x, X)s±(x, A). 
Since the right-hand side of the above equality does not depend on x, we see that 

iy(A)= Cl (l,A) = 0(A). 

Thus the residue of {ss(\ — A) _1 u at A = A„ is equal to 



(3. 



(u(t),s+(t,X n )) C 2dt 



C(x, An) 

0(A n ) 



+ 



(u(t),c(t,X n ))c*dt 



S + (x, An) 
0(A n ) 
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We observe now that the vector-functions s + (x, A n ) and c(x, A„) are collinear, namely, 

s + (-,A n ) = 2 j ' n ^c(-, A n ) = c(-,A„), 
c 2 (l,A n ) y{K) 

and thus expression (I3.8P simplifies to (u, c n )c n /(0(A ri )V'(A n )). Equating this with the 
above expression for the residue of [s^\ — A) _1 u at A = A n , we obtain relation (13.71) for 
a n . The proof is complete. □ 

Proof of Theorem VJ.tA In view of equality (13. 7p and Propositions IA. l l and IA. 2l it suffices 
to prove that the numbers a n := — 0(A n ) and b n := i[)(\ n ) can be represented as 
(— l) n (l + a n ) and (— l) n (l + b n ) respectively, where a n and b n are n-th Fourier coeffi- 
cients of some functions from L p . 

Using formulae (13.31) and (13. 4p . we show that the numbers a n and b n are of the form 

sinA n + J f(s)e iXnS ds 

with f(s) = —isfx(s) G L p (— 1, 1) for a n and f(s) = f2(s) G L p (— 1, 1) for b n . Hence it 
remains to show that (— 1)" sin A n = 1 + e n (gi) for some gi G L p and that 

(-l) n J f(s)e iX » s ds = e n (g 2 ) 

for some g 2 G L p . 

Since by Theorem 13.11 the eigenvalues A n satisfy the relation A n = n(n + |) + A n , 
where A n = e n ((y'i) for some function gx G L p , we find that sinA n = (— l) n cosA n , so 
that 

(-irsinA„-l = ^(-l)^. 

Applying Proposition IA.1I to the element x := (A n ) ne z G X p , we see that there exists 
a function (j^ in L p such that the sum on the right-hand side of the above equality 
equals e n (gi). 

Changing the variables s i— > 1 — It in the integral, we get 

y f(s)e iXnS ds = J 2f(l-2t)e tXn(1 - 2t) dt 

i 

^^ e a„(l-2t) e - 27 rmt dt 

with /(t) := 2i/(l — 2t)e~ mt . Developing the function e* An ( 1_2i ) into the Taylor series 
and then changing summation and integration order (which is allowed in view of the 
absolute convergence of the Taylor series and the integral), we find that 



f(s)e^ s ds = J2lT^(V k f) 
1 k=o ' 



where V is the operator of multiplication by i(l — 2t). In virtue of Proposition IA. II and 
the fact that V has norm 1 in L p , the right-hand side of the above equality gives the 
rz-th Fourier coefficient of some function g 2 from L p , and the proof is complete. □ 
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4. Derivation of the GLM and the Krein equation 

Write c n (-) = c(-,A n ) and a n := ||c„||g 2 . Since the functions {c n } n6 ^ form an 
orthogonal basis of HI, we have 

k 

s-lim V] a n (-,c n )c n = J?, 

n=—k 

where J? is the identity operator in HI. On the other hand, c n = {J* + Jff)v n with 
v n (x) := c (x, A n ) = (cos \ n x, sin Xnx) 1 , so that the previous relation can be rewritten 

as 



(> + JT) [s-lim Y] a n (;v n )v n ]{j? + J{ry = J, 

k~*oo ' * 



n=—k 



which implies 



s-lim V a n (-, v„) Vn = (j + xy\j + JT*) 

n=—k 

Set 

oo 

(4.1) H(s):=V.p. («ne" 2A " sB -e 



-7r(2n+l)s_B 



where for p = 1 the summation is understood in the Cesaro sense, see Lemma 15.21 
Observing that ( • , c (-, A))c (-, A) is an integral operator with kernel 

1 V-X(x-t)B + e -X(x+t)Bj^ 

and that 

k 

s-lim V] (-,v Qri }vo in = 

n=—k 

with v _ n := c (-, Ti (n + |)), we conclude that 

(4.2) (S + JfT)-\S + Jtf*)' 1 = J? + 3?, 

where & is an integral operator in Et, 



(^u)(x) = f F(x,s)u(s)ds 
Jo 



with kernel 

(4.3) F(x,s):=±[H{^) + H(2¥)J]. 

Since J^* is an integral Volterra operator with upper- diagonal kernel, the operator 
J" + X* is invertible and its inverse can be written in the form J? + J^, where X is 
an integral operator in HI with upper-diagonal kernel. By ( 14.2ft one gets 

{J + + &) = {J + JT*)" 1 = J + JT; 

spelling out this equality in terms of the kernels K and F for x > t, we arrive at the 
Gelfand-Levitan-Marchenko (GLM) equation 



(4.4) 



K(x,t) + F(x,t) + / A"(x,s)F(s,t)ds = 0. 
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If Q is continuous, then such is also K, and one has the formula [36j Lemma 12.1.1] 

(4.5) Q(x) = K{x, x)B - BK(x, x) 

relating the potential and the kernel of the corresponding transformation operator. 
This suggests the following algorithm of solution of the inverse spectral problem: given 
the spectral data {(A n ), (a n )}, one constructs first the kernel F via (14. 3 \ and (14.11) . 
then solves the GLM equation (I4.4p for K, and, finally, recovers the potential via (14.51) . 
However, if Q e L p (^#2), then relation ( 14. 5 p becomes meaningless since, by Theo- 
rem 12.31 in this case the kernel K belongs to G p (^2) but can be neither continuous 
nor well defined on subsets of [0, 1] x [0, 1] of Lebesgue measure zero. 

It turns out that for Q G L p (^ 2 ) the restriction R(x,x) of the kernel R of (12. 12|) 
to the diagonal does determine a matrix-function with entries in L p (0, 1). If Q is 
continuous, then fl4.5[) together with (12.121) and the commutator relations 

(4.6) BR(x, t) = R(x, t)B, BJ = -JB 
yields the equality 

(4.7) Q{x) = R(x,x)JB. 

Equation (14. 7p retains sense also for Q 6 L p {^2) as an equality in L p (^ 2 ) and thus 
can be used to recover Q in this situation. 

Our next task is to explain how the kernel R can be determined from the spectral 
data. We do this by deriving below the Krein equation (14.91) . an analogue of the GLM 
equation for R [TUl I3"2"] . 

Applying the operator B to equality (14.41) from both sides, we obtain its counterpart, 

PX 

(4.8) BK(x,t)B + BF(x,t)B + BK{x, s)F{s,t)B ds = 0. 

Jo 

With regard to the commutator relations (14. 6p and BH(x) = H(x)B, we see that 

BK(x, t)B = \ [-R(x, 2=*) + R(x, 2±*) J] , 
BF(x,t)B = ±[-H(¥) + H(2&)] 

and hence 

K(x, t) + BK(x, t)B = R{x, 2±*) J, 
F(x, t) + BF(x, t)B = H(^)J. 

It also follows that 

K(x, s)F{s, t) + BK(x, s)F{s, t)B = \ [R(x, ^)H(s±i)J + R(x, ^)JH{^)} . 

Adding now (14. 4p and (14.81) . combining the above formulae in the resulting expression, 
and using the relation JH(x) = H(—x)J, we arrive at the equation 

R(x,^) + H(^)+ R(x,s)H(^-s)ds = 0, 0<t<x<l. 

Jo 

Subtracting (14.81) from (I4.4p and performing similar transformations, we arrive at the 
above formula with replaced by ^r, and both can now be combined together to 
give 

R(x,t) + H(t)+ I R(x,s)H(t-s)ds = J 0<t<x<l. 
Jo 
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We see that the function R(x,t) '■— R(x,x — t) satisfies the following Krein equation: 

(4.9) R(x,t) + H{x-t) + / R(x,s)H(s-t)ds = 0. 

Jo 

We observe that as soon as a kernel H is given by (14. ip with A n and a n obeying 
the proper asymptotics (guaranteeing that the series for H converges in L p (^# 2 )), the 
integral operator with kernel H(x — t) belongs to the algebra ^,(^#2) introduced in 
Section [2] and J?+J%? is positive in EI (see Lemma l5~4~l) . The Krein equation (14.9P is then 
uniquely soluble and its solution belongs to G p (^2), see Appendix El In particular, 
R(x, 0) = R(x,x) is in L p (0, 1) entrywise indeed. 

Remark 4.1. We notice that the GLM equation (14.41) is the even part of the Krein 
equation (14.91) in the sense that if R is a solution to (j4.9j) . then the function 

(4.10) K{x,t):=i[R{x,^) + R{x,^J)] 



solves (14.41) . Moreover, the condition J? + J$? > implies that the operator J* + & 
is positive in EI and thus, in view of the results of Appendix [Bj guarantees that the 
GLM equation (I4.4p with F of H4.3[) is uniquely soluble for K and the solution belongs 
to G P (J^ 2 )- 

5. Inverse spectral problem 

The purpose of this section is, firstly, to show by limiting arguments that for- 
mula (14. 7p remains valid if the matrix potential Q belongs to L p {^2) and, secondly, to 
justify the algorithm reconstructing the potential Q from the spectral data. Namely, 
we shall prove the following theorem, which constitutes the main result of the paper. 

Theorem 5.1. Assume that {(A n ) ne z, (// n ) n ez} an arbitrary element of SD p; p G 
[l,oo). Then there exists a unique potential Q G L p (^ 2 ) such that (A n ) and (fi n ) 
are eigenvalues of the corresponding Dirac operators =2/1 and respectively. The 
potential Q is equal to R(x,x)JB = R(x,0)JB, where R is the solution of the Krein 
equation (14.91) with H of (14.11) and a n given by (13.71) . 

The reconstruction algorithm proceeds as follows. Given an arbitrary element of SD P , 
we construct functions and ip via relations (I3.6P and then determine the constants a n 
by (13.71) . Since the A„'s and /z n 's interlace, it is easily seen that all a n 's are positive. 
By virtue of the results of [27] there exist functions f\ and f'2 in L p (—1, 1) such that 



(5-1) 



0(A) = cosA + J h(s)e iXs ds, 
if>(\) = sinA + J f 2 {s)e iXs ds. 



Therefore the proof of Theorem 13.31 remains valid, and thus the numbers a n satisfy 
the asymptotics a n = 1 + e n (g) for some g G L p (0, 1). We shall now prove that the 
series (14. ip converges in L p (^ 2 ), so that the function H is well defined and belongs 

to L p (^2). 

Lemma 5.2. Assume that the numbers A n and a n , n G Z, satisfy the asymptotics of 
Theorems \3 . 1\ and \3.3\ for somep G [1, 00). Then the series (14.11) converges in the space 
L p ((— 1, 1), ^#2) (in the Cesaro sense if p — !)■ 
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Proof. Since the matrix B is skew-adjoint and has eigenvalues ±z, it suffices to prove 
that the scalar series 

oo 

(5.2) h{s) := V.p. J2 Ke 2A " iS - e (2n+1) ™) 

n=— oo 

converges in L p (— 1, 1) (in the Cesaro sense if p — 1). We shall justify the convergence 
on (0,1); that on (—1,0) will then follow if we replace a n with and A n with 

— A_i_ n . 

By assumption, a n = 1 + e n (g) for some g G L p , and classical theorems of harmonic 
analysis [2HJ Sec. 1.2] show that the series 

oo 

V.p. J2 K - l)e 2ri7ris 

n=— oo 

converges in L p to g in the required sense. It remains to establish that the series 

oo 

(5.3) V.p. a n {e 2X " is -e {2n+1) ™) 

n=— oo 

is convergent in L p . We shall treat only the case p = 1; the arguments remain the same 
for p > 1 but with partial Cesaro sums replaced by the ordinary partial sums. 

By the definition of summability in the sense of Cesaro we have to show that the 
sequence of partial sums 

TV k 

S N (x):=J2J2 «™(e 2A "" - ^ 2n+l ^ lx ) 

k=l n=—k 

is a Cauchy sequence in L\. Recalling that A n = 7r(n + |) + A„ with A„ = e n (f) for 
some / G L\ and developing e 2Xnlx into Taylor series around (2n + l)nix, we see that 

N k oo /n r . Nm 

[2\ n ix) m 



ml 

k=l n=—k m=l 
oo /„. i„ iV k 



wix (2ix) \ "> \m2mrix 

m=l fe=l n=—k 



oo 



m=l 



v-> (2ia;) , , 
/ — ' m! 



where 



AT fc 



fc=l n=—k 

is the partial Cesaro sum for the function 

/»>:= /*• • •*/ + /*• • •* / *g» m>l 

m m 

(see Appendix [A} . As is well known (cf. [28l Ch. II]), for any / G Li, the partial Cesaro 
sums a N (-,f) converge to / in Li and ||cr JV (-, /)|| ia < \\f\\ Ll - Since 

ll/mlk^a + lblMII/IISi, 
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the Lebesgue dominated convergence theorem shows that the limit 



CO 



lim S N (x) = e-™ V ^-/ m (x) 

m=l 

exists in L±, whence the series (15.31) converges in L\ in the Cesaro sense. The lemma 
is proved. □ 

Next we show that the GLM equation (14.41) and the Krein equation (14. 9 j) have unique 
solutions that belong to the space G p (^# 2 )- To this end it suffices to show (see Appen- 
dix [B] and Remark 14.11 for details) that the operator J? + M 3 is positive in H, where 
3ft 3 is a Wiener-Hopf operator given by 

JTu(x) := [ H(x - s)u(s) ds. 
Jo 

As a preliminary, we show that certain systems of functions form Riesz bases in L 2 
or EL 

Lemma 5.3. Assume that the sequence (A n ) is as in Theorem \5.1\ Then 

(a) the system {e lXnX ) n & forms a Riesz basis of L 2 {— 1, 1); 

(b) the system (v n ) ne g ; with v n := c (-, A n ) ; forms a Riesz basis o/EL 

Proof. In view of relation (13.31) the numbers A^ are zeros of the exponential function of 
sine type with indicator diagram [—1,1] |34j . so that item (a) follows from [6j Propo- 
sition II.4.3]. 

Assume now that g := (g\, t7 2 )* is an arbitrary element of EL We extend g\ and 
g 2 to functions gi and g 2 on (—1, 1) in the even and odd way, respectively, and set 
g := gi + ig 2 . By (a), there exists a unique sequence (a n ) ng z in £ 2 , for which 



~g = a ^ tX " X 



n=—oo 



the series being convergent in L 2 (— 1, 1). Taking the even and odd parts of the above 
equality, we arrive at the relations 



CO 



9i( x ) = a "Cos(A„x), ig 2 {x) = i ^2 a n sm(X n x), 

n=—oo n=— oo 

i.e., at the representation 

oo 

g = ^ a „v„ 

n=— oo 

in EL We observe that there is a constant C independent of g such that 

oo oo 

c- 1 J2 Kl 2 <ll^ll! 2 = 2||g||^<c J2 K\ 2 - 

n=— oo n=— oo 

Hence the system (v n )„ 6 ^ is a Riesz basis of EL □ 

Lemma 5.4. Assume that {(A n ), (/i n )} an arbitrary element of SD P and inai t/te 
function H is constructed as explained above. Then the corresponding operator J? + 
is positive in EI. 
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Proof. We notice that the matrix B is skew-adjoint and has eigenvalues ±z. It is clear 
that the subspaces ker(I? ± i) eg) L 2 (0, 1) are invariant subspaces of Jff '. Therefore ffl 
is unitarily equivalent to the direct sum © J^L, where 

(M±g)(s)= h(±(x - s))g(s)ds, g e L 2 , 
Jo 

and h is the function of (15.21) . Thus positivity of J' + is equivalent to that of both 
/ + J£+ and / + M'- . We shall prove only that the operator / + M'^ is positive in L 2 ; 
the positivity of the other one is established analogously. 

Since the systems ( e ^ 2n+1 ^ 7 " s ) ngZ an d ( e2A "* s ) ng x constitute respectively an orthonor- 
mal and a Riesz basis of L 2 (cf. Lemma [5.31 and P, Sect. II. 4. 2]), it is easy to see that 

k 

((J + J^)/,/) = (/,/)+lim [«nl(/ ; e 2A "-)| 2 -|(/,e^ +1 ^)| 2 ] 



n=— k 



2A n is\ 1 2 



= <* n \(f,e 

n=—oo 

Since the numbers a n are uniformly bounded away from zero, we conclude that there 
is a C > such that ((/ + Jf + )f, f) > C\\f\\ 2 for all / £ L 2 , i.e., the operator / + 
is (uniformly) positive in L 2 . □ 

According to Theorem IB.21 positivity of the operator J? + M 3 implies that the 
GLM equation (14.41) with F given by (14.31) and the Krein equation (14.91) have unique 
solutions K,R £ G p {^£ 2 ). We have shown in Section 0] that in the smooth case the 
solution K is the kernel of the transformation operator J* + Jff for the pair (&/q, 
with Q := R(-, 0)JB. Based on this result, we treat here the general case by a limiting 
procedure. 

Theorem 5.5. Assume that {(A n ), (fi n )} is an arbitrary element of SD P and that H £ 
L p {^ 2 ) is a function of (14.11) constructed as explained above. Let also K and R be the 
solutions of the GLM equation (14.41) and the Krein equation (14. 9p respectively. Denote 
by the integral operator with kernel K. Then there exist a unique Q £ L p (^ 2 ) — 
namely, Q = R(-, 0)JB — such that the operator J? + is the transformation operator 
for the pair s^q and &/o- 

Proof. We shall approximate the function H in the norm of L p {^ 2 ) by a sequence 
(-f^)^i °f real-valued, smooth (say, infinitely differentiable) ^# 2 -valued functions so 
that the following holds: 

(a) for every / £ N, the GLM equation (I4.4p with H replaced by Hi has a unique 
solution Ki, and the corresponding integral operators J£f converge to as 
I — > oo in the uniform operator topology of H; 

(b) for every I £ N, there exists Qi £ L p {^ 2 ) of the form (II. ip such that J? + 
is a transformation operator for the pair s^q 1 and j^ Q ; 

(c) the matrix-functions Qi converge to Q := R(-,0)JB in L p {^ 2 ). 

If (a)-(c) hold, then by Theorems 12 . 3 1 and 12 . 41 1 he operators J? + converge in %(^# 2 ) 
(and hence in the uniform operator topology of H) to an operator J? + J^q, which is 
the transformation operator for the pair (jz/q, sz?q). Thus ,Jff = J£q yielding the result. 
The uniqueness of Q is obvious. 
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The details are as follows. Using (A n ) and (/i n ), we construct the sequence of con- 
stants a n and set 

i 

Hl (s) := (ane~ 2X " sB - e -^ 2n+1 ^) 

n=—l 

(for p — 1, we replace Hi by the corresponding partial Cesaro sum | Ylk=i ^fc)i then 
Hi ^ H in L p {^2) as I — > oo by Lemma [5. 2[ 

We observe that this choice of if/ corresponds to setting a n = 1 and A„ = 7r(n + |) 
for all n with |n| > I, so that in view of Lemma 15.41 the Wiener-Hopf operators Jtf 
with symbol Hi satisfy the condition J? + M{ > 0. Hence by Corollary IB.3I and 
Remark 14. II the GLM equation (14. 4p with H replaced by Hi has a unique solution Ki. 
This solution Ki belongs to G p (^ 2 ), and hence the corresponding integral operator 
is bounded. Since the relation Hi — > H in L p {^2) as / — > oo implies that J£? — > J^ 7 
in ^,(^# 2 )) we conclude that — > as / — > 00 in the topology of the space ^(^# 2 ) — 
and thus in the uniform operator topology of H, see Appendix [B] This establishes (a). 

Moreover, by the well-known result for continuous potentials |36| Ch. 12.3-4] the 
operator + is the transformation operator for the pair {&?q 1 , £/o) with Qi(x) := 
Ki(x,x)B — BKi(x,x). As was explained in Section HJ this yields the relation Qi(-) = 
Ri(-, 0)JB, where Ri is the solution to the Krein equation (14. 9 p for H = Hi. Thus (b) 
is fulfilled. 

To prove (c), we observe that, according to the results of Appendix [Bj the solution 
R to the Krein equation (14.91) depends continuously in the norm of the space G p {^2) 
on the function H e L p (^ 2 )- Therefore the sequence i?;(-,0) converges in L p (^# 2 ) to 
the function R(-,0). The proof is complete. □ 

The last step of the reconstruction procedure is to show that the numbers X n and 
fi n we have started with are the very eigenvalues of the operators srf\ and with the 
potential Q just found. 

Since the solution K to the GLM equation (14. 4p generates a transformation operator 
.y+J^ for the pair (^q, ^0), the functions c(-, A) := (J?+,J^)c (-, A) belong to domi^Q 
and satisfy the relation £/qc(-, A) = Ac(-, A). We set := c(-, X k ) and show that these 
functions are orthogonal and that the a k are the corresponding norming constants. 

Lemma 5.6. The system of functions {ck}kez is an orthogonal basis ofW. Moreover, 
for the above numbers a k (defined at the beginning of this section), we have 

(c fc ,Q) = a^5ki, 

where 5 k i is the Kronecker delta. 

Proof. Denoting by v n the function Co(-, A n ) and recalling that the integral operator & 
with kernel F of g3J is related to by {J + + + JT)* = J (see 

Appendix IBl for details), we conclude that 

(c k , ci) = ((y + + JT)v fc , v,) = ((S + ^)- x v fe , v,). 

Reverting the arguments of Section HJ we see that the operator J? + & is equal to 

k 

s-lim ) a n (-,v n )v n . 

n=—k 
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Since the sequence (v n ) ng z is a Riesz basis of H in view of Lemma 15.31 the inverse of 
J? + & can be represented as 



n 



(J? + 3?) ^s-lim V a m 1 (-,v m )v m , 

n— >oo ' ' 

m=—n 

where (v m ) is a basis biorthogonal to (v m ) (see [HJ Ch. VI]). Therefore 

n co 

{{J + &)~ l \ k , vi) = s-lim a~ 1 (vji ! ,v m )(v m ,v i ) = } oc^8 km 8 m i = a^Ski- 

n. — >nn t * f * 



n— >oo 

m= 



Completeness of the system {c n } nG z immediately follows from the fact that the system 
{vnjngz is complete and that + J?T is a homeomorphism of H, and the lemma is 
proved. □ 

Next we show that the numbers X n are indeed the eigenvalues of the operator s$\. 
According to what was said above, it suffices to show that ci(l,Afc) = 0. For the 
operator &/q, one has the Green formula 

(**qc{; A),c(-,/z)) - ((c(-, \),£/ Q c(-,n)) = c 2 (l, A)ci(l,/x) - ci(l, A)c 2 (l, //); 

taking therein A = A& and /i = A n and using the previous lemma, we arrive at the 
relation 

(5.4) c 2 (l, A fc )ci(l, A n ) = ci(l, A fe )c 2 (l, A n ). 

Assume first that none of the numbers ci(l, A&) vanishes. Relation (15.41) implies that 
there is a constant 7 such that, for all k G Z, 

c 2 (l, A fc )/ci(l, A fc ) = 7. 

Then are eigenvectors corresponding to the eigenvalues A& of the operator that is 
the restriction of s^q by the boundary condition w 2 (l) = 71*1(1). In other words, the 
numbers A n are zeros of the function 

c 2 (l, A) - 7Ci(l, A) = sin A - 7 cos A + J f(s)e lXs ds 

for some / G L p (—1, 1). However, the standard arguments based on Rouche's theorem 
(see, e.g., [HI Ch. 1.3]) show that the zeros \ n of the function c 2 (l, A) — 7Ci(l, A) obey 
the different asymptotics X n = nn + arctan7 + o(l), which leads to a contradiction. 

Therefore there is a k G Z such that ci(l,Afc) = 0. Then c 2 (l,A&) 7^ due to the 
uniqueness of solutions to the equation £q(vl) = A^u, and relation (15.4)) implies that 
c 2 (l,A„) = for all n G Z. In other words, the functions c n are the eigenvectors of 
the operator srf\ corresponding to the eigenvalues A n . Since by Lemma [5.61 the system 
{c n } n( zz is complete in HI, the operator srf\ has no other eigenvalues. 

It remains to prove that \x n are the eigenvalues of s&2- We denote by ft n these eigenval- 
ues and construct the corresponding function ijj of (I3.6p . Recalling expression (13.71) . we 
conclude that ip(\k) = V'(Ajfc) for all k G Z. Since the function ip has the representation 



$(A) = sinA + J f 2 {s)e lXs ds 
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for some f 2 G L p (— 1, 1) (see the proof of Theorem 13.11) and ip has a similar repre- 
sentation with some / 2 G L p (— 1, 1) instead of / by (15. ip . we see that the function 
lu := f'2 — /2 is such that 



for all n G Z. Recalling that the system of functions {e* AnS }„ e z is closed L p (— 1, 1) 
(this follows from [351 Theorem III] for p > 1 and from [23, Lemma 3.3] for p = 1), 
we conclude that lu = 0. Thus the numbers u n = Jx n are the eigenvalues of the opera- 
tor s&ii and the reconstruction procedure is complete. As the spectral data determine 
the function H (and thus the transformation operator J? + Jtf) unambiguously, the 
potential Q is unique. The proof of Theorem 15. II is complete. 

We observe that, in passing, we have solved the inverse spectral problem of re- 
constructing the potential Q of the Dirac operator from the spectrum of stf\ and the 
corresponding sequence of norming constants. Namely, the following is true. 

Theorem 5.7. Sequences of real numbers (X n ) n( z^ and positive numbers (a n ) nG z are 
respectively the sequences of eigenvalues and norming constants of an operator for 
some Q G L p {^ 2 ), V G [1, 00), if and only if the following holds: 

(i) the numbers \ n strictly increase and obey the asymptotics of Theorem \3.1\ 

(ii) the numbers a n > obey the asymptotics of Theorem \3.3[ 

If (i) and (ii) hold, then Q is given by R(x, x)JB = R(x, 0)JB, where R is the solution 
of the Krein equation (14.91) with H of (14.11) . 

The reconstruction algorithm proceeds as the previous one, except that the first and 
the last step (related to the spectrum (// n )) should be omitted. 

In a similar manner we can also treat the inverse spectral problem for the operator stf 2 
(or operators corresponding to arbitrary separated boundary conditions). 
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Appendix A. Fourier transform in L p (0, 1) 
For any / G L p (0, 1), we denote by e n (f), n G Z, its n-th Fourier coefficients, i.e., 



The vector space X p is algebraically embedded into ^oo(Z) and becomes a Banach space 
under the induced norm 



For any x = (x n ) and y = (y n ) in (Z) we shall denote by xy the entry wise product 
of x and y, i.e., the element of ^oo(Z) with the n-th entry x n y n . 




uj(s)e iXnS ds = 





e(f)\\x p := 
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Proposition A.l. X p is a commutative Banach algebra under the entrywise multipli- 
cation, i.e., 

(A- 1 ) ll x y|U P < ||x|UJy|U P - 

Indeed, inequality ( 1A.1I) follows from the fact that e n (f)e n (g) = e n (f * g), where 

if*9)(x) ■= [ f(x-t)g(t)dt 
Jo 

is the convolution of / and g (f being periodically extended to (—1, 1) by f(x + 1) = 
f{x), x G (0, 1)), and from the inequality 

\\f *q\\l p < H/IUJMIv 
The following statement is an analogue of the well-known Wiener lemma. 

Proposition A. 2. Assume that f G L p (0, 1), where p G [1, oo). If 1 + e n (f) ^ for 
all n6N, then there exists a function g G L p (0, 1) such that 

(l + enif))' 1 = l+e n (g). 

Proof. To begin with, we adjoin to X p the unit element 5 (with all components equal 
to 1) and denote the resulting unital algebra by X p . Assume that the assumptions of 
the lemma hold and denote by x an element of X p with components x n := 1 + e n (f). 
We shall prove below that x is invertible in X p ; since e n (f) —>■ as n — > oo, this will 
imply that x^ 1 = 5 + y for some y G X p as required. 

As is well known [16] , the element x is invertible in the unital Banach algebra X p if 
and only if x does not belong to any maximal ideal of X p . Proceeding by contradiction, 
assume that there exists a maximal ideal m of X p containing x. Since X p includes all 
finite sequences and none of x n vanishes, m also contains all finite sequences. Finite 
sequences form a dense subset of X p because the set of all trigonometric polynomials is 
dense in L p (0, 1). Recalling that maximal ideals are closed, we conclude that X p C m. 
Next we observe that X p is a proper subset of m (e.g., x belongs to m \ X p ) and that 
X p has codimension 1 in X p . Hence m = X p , which contradicts our assumption that 
m is a maximal ideal of X p . As a result, x is not contained in any maximal ideal of X p 
and thus is invertible in X p indeed. The lemma is proved. □ 

Appendix B. The GLM equation and factorisation of Fredholm 

OPERATORS 

In this appendix, we shall explain relationships between solubility of the GLM equa- 
tion and factorisation of related Fredholm operators in some special algebras. We refer 
the reader to the books [HI [20] for related concepts and basic facts. 

Write H := L 2 ((0, 1),C 2 ) and denote by 9S (by SS^) the Banach algebra of all 
bounded (compact) operators in H. Denote also by P t , t G [0, 1], the operator in HI of 
multiplication by X[o,t], the characteristic function of the interval [0, t). Set 

@to ■= {B e -SSoo | Vt G [0, 1], P t B(I - P t ) = 0}, 
:= {B G | Vt G [0, 1], (/ - P t )BP t = 0}; 

then are closed subspaces of and ^+ H = {0}. We also observe that the 
operators in SS^ are Volterra ones. 
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Recall that ^(-#2), P G [1, 00), stands for the algebra in of integral operators 
over (0, 1) with kernels in the space G v {y$-i) introduced in Section [2J The sets 

% + (^2) = %{^2) n %{Jt*) ■= %{^2) n ^" , 

are subalgebras of ^,(^# 2 ) consisting of operators with lower- and upper-triangular 
kernels respectively and £f p (^# 2 ) = % + (^) + ^ p ~(^ 2 )- 

We say that an operator + Jzf, Jzf G £^(^2), admits a factorization in ^(^#2) if 

(B.i) j + ^ = + + ^-y 1 

with some G Sfjf 

The following two theorems were established in [131 HI] for the space L 2 (0, 1), how- 
ever, their generalisation to our situation is straightforward. 

Theorem B.I. If J? + Jzf admits a factorization in ^p(«/# 2 ), then the operators J^ A± = 
J^ ± (^f) are unique. Moreover, the set $ p of operators Jzf G ^,(^2), /or which the 
operator J? + Jzf is factorisable, is open in %(y# 2 ), functions 

are continuous. 

Theorem B.2. Assume that Jzf G ^,(^# 2 ). For i/ie operator J? + Jzf to admit a 
factorisation in ^(^2), it necessary and sufficient that the operators I + P t J£Pt 
/iai>e a trivial kernel in HI /or eac/i t G [0, 1]. 

We remark that for a self-adjoint operator Jzf G ^,(^2) the requirement that the 
operators I+Pt J^Pt have a trivial kernel in H for all t G [0, 1] is equivalent to positivity 
of + Jzf in H. 

Assume that ^+J£ is factorisable in ^ p (^# 2 ), so that (IB.ll) holds. Applying 
to both sides of this equality and using the fact that (J^ + J^ - ) -1 = J? + for some 
G =5$"", we derive the relation 

(B.2) jr + + ^ + j^ + ^ + pr+^) =0, 

where ^ + denotes the projection operator of 5^,(^# 2 ) onto £f p + (^# 2 ) parallel to ^,~(^# 2 ). 
This relation is an abstract analogue of the Gelfand-Levitan-Marchenko (GLM) equa- 
tion; indeed, in terms of the kernels K + and L of the operators + and Jzf we get 

(B.3) K + (x,t) + L(x,t) +/ K + (x,s)L(s,t)c/s = 0, 0<t 
cf. (TOD and (jUD - 

We see that if an operator J? + Jzf is factorisable in ^,(^#2), then the abstract 
GLM equation flR2l) has a solution JT+ = JT + (j5f) G Sf+(-# 2 ). Conversely, if JT+ G 
£f+(^f 2 ) is a solution of equation flR2j) . then ^ + (J^+ + Jzf + JT+Jz^) = 0, i.e., X := 
JT+ + Jzf + belongs to 9f(jtr 2 ), and the relation 

{J + JT+)(^ + Jz^) = J + 

holds. Since J + X has the form + for JT~ := (^ + ^0^ - ^ G Sf-(-#f 2 ), 

we conclude that J? + Jzf is factorisable in ^ p (^# 2 ). 

Summing up, we derive the following assertion on the solubility of the GLM equa- 
tion (El . 



< x < 1, 
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Corollary B.3. Assume that ££ is a selfadjoint operator in ^,(^#2) with kernel L such 
that J* + S£ is positive. Then equation (IB.3I) is uniquely soluble, and the solution K + 
belongs to (j p (^# 2 ) an d depends continuously therein on L 6 G p (^# 2 )- 
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